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Abstract 

We investigate both theoretically and experimentally the rotational dynamics of micrometric 
droplets of dye-doped and pure liquid crystal induced by circularly and elliptically polarized laser 
light. The droplets are dispersed in water and trapped in the focus of the laser beam. Since 
the optical torque acting on the molecular director is known to be strongly enhanced in light- 
absorbing dye-doped materials, the question arises whether a similar enhancement takes place also 
for the overall optical torque acting on the whole droplets. We searched for such enhancement 
by measuring and comparing the rotation speed of dye-doped droplets induced by a laser beam 
having a wavelength either inside or outside the dye absorption band, and also comparing it 
with the rotation of pure liquid crystal droplets. No enhancement was found, confirming that 
photoinduced dye effects are only associated with an internal exchange of angular momentum 
between orientational and translational degrees of freedom of matter. Our result provides also the 
first direct experimental proof of the existence of a photoinduced stress tensor in the illuminated 
dye-doped liquid crystal. Finally, peculiar photoinduced dynamical effects are predicted to occur 
in droplets in which the molecular director is not rigidly locked to the flow, but so far they could 
not be observed. 
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I. INTRODUCTION 



It has been known since the early eighties that hght can transfer its angular momentum to 
liquid crystals with high efficiency, causing the rotation of the average molecular orientation, 
as specified by the molecular director [1]. 

In the early 90's it was discovered that adding small amounts of certain dyes to the liquid 
crystals, the light-induced torque on the molecular director could be greatly enhanced |2l [3] . 
This effect, initially rather puzzling, was later explained by a model relying on the reversible 
changes of intermolecular forces occurring between photoexcited dye molecules and liquid 
crystal host [4J. This model was subsequently extended and refined [S] and con&med with 
several independent experiments |6l [71 El El flO] - 

The understanding achieved with this model has also provided an answer to a fundamental 
question related to the observed torque enhancement: where does the additional angular 
momentum come from? Clearly, the angular momentum that is transported by the optical 
wave impinging on the material cannot be affected by the presence of dye. It is true that 
absorption of light due to the dye does lead to some additional transfer of angular momentum 
from light to matter and therefore to a variation of the light-induced torque. However, this 
effect cannot account for the magnitude of the observed torque enhancement [6] and for 
its peculiar dye-structure dependence [7J. The clearest solution to this puzzle was put 
forward in Refs. [5] and [TTJ and is as follows. The angular momentum transfer from light 
is indeed approximately unchanged. Light absorption however triggers a transfer of angular 
momentum between different internal degrees of freedom of the liquid crystal, namely from 
the center-of-mass molecular degrees of freedom, corresponding macroscopically to fluid flow, 
to the molecular-orientation degrees of freedom, corresponding to the molecular director. 
Being an internal transfer, the total angular momentum given to the material as a whole is 
not affected by the presence of dye. Another consequence of this internal transfer is that, 
by the force action-reaction principle, in dye-doped liquid crystals there should be another 
(opposite) torque acting on the fluid in the form of a photoinduced stress-tensor [5j. 

However plausible they may be, both these predictions have never been tested directly in 
an experiment. Here we address these questions by studying the rotational behavior of pure 
and dye-doped liquid crystal droplets of micrometric size, optically trapped in water, under 
the effect of circularly and elliptically polarized infrared and visible laser light. In particular. 
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we measured the rotation speed of dye-doped droplets illuminated with visible light, having 
a wavelength in the dye absorption domain, and compared it with that of undoped droplets 
and with the case in which the illumination is by infrared light, not absorbed by the dye. 
Should we observe some significant dye-induced enhancement (or other anomalies) of the 
droplet rotation speed, this would imply that the current understanding of the photoinduced 
torque as entirely due to an internal exchange of angular momentum is not correct, or at 
least not complete. 

But what do we exactly mean by droplet rotation? Prom a theoretical point of view, a 
liquid crystal droplet actually has two independent rotational degrees of freedom: the average 
molecular director, corresponding to the average orientation of the liquid crystal molecules 
in the droplet, and the internal rotational flow, corresponding to a fluid-like motion of the 
molecule centers of mass (this is already a strong simplification: a liquid crystal droplet 
has actually an infinite number of rotational degrees of freedom associated with changing 
inhomogeneous configurations of director and fiow) . These two rotational degrees of freedom 
of the droplet are always strongly interacting with each other, via viscosity, elastic forces, 
and possibly also photoinduced effects. This interaction may be strong enough to lock the 
two degrees of freedom together into a single one, so that the droplet effectively behaves 
as a rigid body. This is what we actually observed in our experiments. Nonetheless, by 
properly choosing the experimental parameters, it should be possible to reach a regime in 
which the two rotational degrees of freedom become effectively unlocked. We would expect 
that in this "unlocked" regime, photoinduced dye effects become important, even if the 
angular momentum exchange is entirely internal. We investigated theoretically this case 
and predicted the occurrence of highly non-trivial photoinduced rotational effects, which we 
then tried to observe experimentally. 

This article is organized as follows. The general continuum theory of the light-induced 
dynamics of dye-doped liquid crystals is reviewed in Sec. II. In Sec. Ill, we then apply 
this theory to the case of a spherical droplet immersed in water. With the help of several 
strongly-simplifying approximations, analytical solutions are found, predicting what should 
be the light-induced rotational behavior of the droplets as a function of light power and 
polarization and of droplet size. Our experiments are then described and discussed in Sec. 
IV. Our results are finally summarized in Sec. V. 
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II. GENERAL THEORY OF PHOTOINDUCED DYNAMICS OF NEMATICS 



A. Dynamical fields and equations 

In general, the dynamics of a nematic liquid crystal under the action of laser light is 
defined by the temporal and spatial dependence of the following fields: (i) the molecular 
director n specifying the local average molecular alignment [12], (ii) the velocity vector v 
defining the flow of matter, and (iii) the electric and magnetic fields of the optical wave, E 
and B, respectively. 

The optical fields E and B are governed by the usual electromagnetic Maxwell equations 
in (anisotropic) dielectric media [13]. The two material fields n and v are respectively 
governed by the director torque balance equation. 



and by the Newton equation for the acceleration (or momentum conservation law). 



where dj stands for the partial derivative d/dxj and the usual sum convention over repeated 
indices is understood. In these equations, t*°* is the total torque density acting on the 
molecular director, t*°* is the total stress tensor associated with a fluid displacement with 
no director rotation (this corresponds to Ericksen's definition of stress tensors in liquid 
crystals [14J), J„ is a moment of inertia per unit volume associated with the rotation of 
the nematic director, which is actually negligible in all practical cases (it is included only 
for making the equation physical meaning clearer), and p is the mass density, hereafter 
assumed to be constant (incompressible fluid approximation). Moreover, in both equations 
above we have used the so-called "material" or "convective" time derivative, defined as 
d/dt = d/dt + Vjdj, corresponding to a derivative taken following the fluid element along its 
motion. Incompressibility actually adds the following constraint on the velocity field v: 




(1) 



P 



dvi 
~dt 




(2) 



diVi = 0, 



(3) 



which is a particular case of the mass continuity equation. 
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B. Constitutive equations 



Equations ([T]) and (|2]) are to be completed with the appropriate constitutive equations for 
the total torque density and stress tensor. To first order in all deviations from equilibrium 
and in all gradients, we may distinguish five additive and independent contributions to both 
the torque density and the stress-tensor [m [15]: hydrostatic pressure (hp), elastic (el), 
viscous (vis), electromagnetic (em), and photoinduced (ph), the latter being that associated 
with dye effects. In summary, one may write 

r*°t = ^r" (4) 

and 

a 

Each of the ten t° and {a = hp, el, vis, em, ph) terms appearing in these two equations 
has a well defined (first-order) constitutive dependence on the dynamical fields, which we 
will now briefiy discuss. 

First, there is actually no torque density associated to pressure effects, i.e., t^^ = 
identically. The hp stress tensor term has instead the usual simple form t^f = —pSji, where p 
is the pressure field and 6ij is the Kronecker delta. In the incompressible-fluid approximation 
we are adopting, p must be treated as a pure "constraint force", i.e., assuming just the space 
and time dependence needed to ensure continuous validity of Eq. 

Next, the constitutive laws of the elastic and viscous torque densities and stress-tensors 
are fully standard and we refer to Refs. [H] or (TS] for their explicit expressions. 

Let us now turn to the electromagnetic terms. They are also standard, but it is nev- 
ertheless convenient to introduce them explicitly here. We assume for the time being that 
the electromagnetic flelds present in our systems are associated only with an approximately 
monochromatic optical wave having a given vacuum wavelength A. Let us first introduce 
the uniaxial optical dielectric tensor of the liquid crystal 

Eij = eo{e±6ij + eaUinj), (6) 

where Eq is the vacuum dielectric constant, 6± the relative dielectric constant for E ± n and 
Ea the relative dielectric anisotropy. Neglecting all magnetic effects at optical frequencies, 
the electromagnetic torque density is then given by 

T^- = (DxE) = £o£,((n-E)(nxE)), (7) 
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where D is the usual dielectric displacement field {Di = SijEj) and where () denotes a 
time average over an optical cycle. The electromagnetic stress tensor can be written in the 
following form [131 [IS] : 

where /xq is the vacuum magnetic permittivity constant and Fem is the electromagnetic free 
energy density at given electric field E and magnetic field strength H = B//io, which in our 
case can be written as 



~ _ 1 52 



2 2/io 



£0 B'^ 
= -— {e±6ij + eaUinj) EiEj - - — . (9) 
z Zfio 

Let us finally consider the photoinduced terms, i.e., those appearing in dye-doped liquid 

crystals when illuminated with light having a wavelength falling within the dye absorption 

band. These are absolutely nonstandard. In the limit of small light intensities, simple 

symmetry arguments show that the photoinduced torque density must be identical to the 

electromagnetic one except for the replacement of the dielectric anisotropy ea with a new 

material constant (, proportional to the absorbance (or to dye concentration) [5]. Therefore, 

its explicit expression can be written as follows: 

TP'^ = Ceo((n-E)(nxE)) (10) 

(the constant Eq is inserted for making ( dimensionless). Similarly, we may apply symmetry 
arguments for identifying the most general possible expression of the photoinduced stress 
tensor. This expression contains seven new unknown material constants (all proportional to 
the absorbance): 

tf=(^[aiE^ +a2(n-E)2] S,j +a^EiEj + [a^E^ 

+a5(n-E)2 niUj +{u.-'E){aQEinj +a'jEjni)y (11) 

The effects of this stress tensor have never been measured or even detected in ordinary 
liquid crystals, although related photoinduced effects may have been observed in polymeric 
nematic elastomers [T71 118] . 

We have thus completed the set of constitutive equations needed to close the dynamical 
equations ([T]) and (|2]). 



C. Angular momentum conservation 



Before moving on to the specific case of a droplet of liquid crystal immersed in water, it 
is convenient to see how the law of angular momentum conservation enters our problem. 

In contrast to the case of the conservation of (linear) momentum, which provides an 
additional dynamical equation, the law of angular momentum conservation actually sets 
only a general constraint on the possible constitutive laws of torque densities and stress 
tensors. For any given volume V of material, the corresponding angular momentum rate of 
change will be given by the following law: 

^ = M*°\ (12) 
where L is the total angular momentum within volume V and M*°* is the total external torque 



acting on it. The two sides of Eq. (12) can be deduced by multiplying the corresponding 
sides of Eq. ^ vectorially by r, integrating them over the volume V, and then adding to 
them the volume integral of the corresponding sides of Eq. ([T]). In this way we obtain 

L = l^iprxY + Innx^)dV (13) 

and M*°* = Ea M", with 

Mt= [ ie,,,xjdktt, + rndV, (14) 

J V 



where eijh is the fully antisymmetric Levi-Civita tensor [19]. However, in order for Eq. (12) 
to be equivalent to a local conservation (continuity) law, it should be possible to reduce 
the total external torque to a pure surface integral over the boundary dV of V, such as the 
following: 

Mt = i ei,H{xjtt, + nj5t,)dA,, (15) 

•JoV 

where s^^ is a new material tensor expressing the torque per unit area exchanged by the 
director n directly through the surface and dAj. denotes a vector having direction equal to 
the local surface normal (pointing outward) and modulus equal to the area of the surface 
element [TitfT^]. 



By equating the two expressions (14) and (15) of the external torque for any possible 



volume V, and exploiting the standard divergence theorems, one obtains the local identity 



i - eijht^h + dkieijhTijS'^f,). (16) 
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Generally speaking, this identity holds true only for a = tot, and not separately for each term 
a = hp, el, vis, em, ph. However, since in our first-order theory these five terms can be tuned 



independently from each other, identity (16) must hold true also for a = hp, el, vis, em, ph, 
separately. 

Moreover, since in the first order approximation only the elastic forces (torque density 
and stress-tensor) are taken to depend on the director spatial gradients, we may deduce 



from Eq. (16) that s^^ is nonzero only for the elastic contribution a = el. In all other 
cases one must have s" = within first order approximation. Therefore in these cases, 
owing to angular momentum conservation, the stress tensor defines completely the torque 
density, or conversely, the torque density defines the antisymmetric part of the stress tensor. 



In particular, Eq. (16) with a = ph yields the following relationship between the material 



constants appearing in expressions (10) and (11): 



(eo = aj - ttQ. (17) 

In concluding this Section, we note that the theory we have just described and in par- 
ticular the constitutive equations we have adopted for the photoinduced torque density and 
stress tensor already imply that the photoinduced angular momentum transfer associated 
with the dye is fully internal to the liquid crystal. Indeed, the flux of angular momentum 



through any boundary surface, as given by Eq. (15), will vanish identically if the stress ten- 
sor vanishes on it. Since we have assumed that all material constants a, are proportional to 
the dye concentration, they should vanish on a surface lying just outside the liquid crystal 
(in water) and therefore the flux of angular momentum through such a surface will vanish, 
i.e., one has M^^ = 0. We note that the same argument does not hold for M'^™, as f^™ has 
a finite value also in the isotropic liquid. 

Is there a possible way out of this conclusion, justifying perhaps a hypothetical photoin- 
duced flow of angular momentum out of the liquid crystal? Within a first order theory of 
the constitutive equations the answer is no. However, it cannot be a priori excluded that 
higher-order terms in the constitutive equations become important in specific situations and 
justify a strong exchange of angular momentum with the outside. For example, at the sur- 
faces between the liquid crystal and the surrounding medium, the mass and composition 
densities suffer sharp discontinuities. Therefore, a first-order theory in the spatial gradients 
is clearly not justified anymore (this, by the way, is just how surface anchoring enters the 
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problem). It would then be conceivable that higher-order terms in the photoinduced torque 
density and stress-tensor expressions could give rise to photoinduced surface effects leading 
to a significant angular momentum exchange with the outside, i.e., to a NP^ ^ 0. In the 
end, this hypothesis can only be tested, and eventually ruled out, experimentally. 

III. DROPLET ROTATIONAL DYNAMICS 

Equations ([T]), (|2]), and Eq. (|3]), supplemented with all constitutive equations for torque 
densities and stress tensors, completely define the light-induced dynamics of the liquid crys- 
tal. In the case of a droplet of liquid crystal immersed in water one should also include in 
the system the appropriate boundary conditions at the droplet surface. We limit ourselves 
to mentioning them: continuity of fluid velocity and forces across the boundary, continuity 
of tangential components of E and of normal components of D, continuity of B, and ap- 
propriate anchoring conditions on n. Moreover, one should account for the dynamics of the 
water surrounding the droplet. The latter is also governed by Eq. ([2]), but with a simpler 
expression of the stress-tensor, including only hydrostatic pressure, Newtonian viscosity and 
the electromagnetic stress tensor. 

The resulting system of equations is clearly very complex and an exact solution can be 
determined only numerically. In the following, we instead approach the problem analytically 
with the help of several approximations. 

First, we will assume that the droplet is always perfectly spherical, with a radius R and 
a total mass m [spherical droplet approximation, SDA). This is an approximation because 
anchoring effects combined with elastic interactions may actually slightly distort the shape 
of the droplet into an ellipsoid. The molecular director configuration within the droplet will 
be taken to be axial, namely the director n has a well defined uniform direction close to the 
center of the droplet, while it will be distorted to some extent close to the surface due to 
anchoring [20t |21] |2^ . This allows one to properly define an average molecular director of 
the droplet. 

Second, for calculating exactly the overall electromagnetic torque acting on the droplet 
we would have to determine in all details the light propagation within the spherical droplet, 
including all the birefringence and wave diffraction (Mie scattering) effects: an exceedingly 
complex task. Instead, following a common practice in the literature [221 |23l [21], we will use 
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an approximate expression of the electromagnetic torque obtained by simply replacing the 
spherical droplet with a homogeneous slab of liquid crystal having a thickness equal to the 
droplet diameter and the strongly focused light beam with a plane wave {planar symmetry 
approximation, PSA). This approximation will tend to become more exact in the limit of 
large droplets and weakly-focused light beams. 

Third, we will restrict the possible dynamics of the fluid and the director to either one 
of the following two approximate models: (i) the droplet behaves exactly as a rigid body, 
i.e., rotating only as a whole and with the director perfectly locked to the fluid {rigid body 
approximation, RBA); or (ii) the droplet fluid flow and director are allowed to have different, 
although uniform, rotation dynamics but the director field is taken to be perfectly uniform 
{uniform director approximation, UDA). 

Let us now go into the details of the outlined approximations. The SDA approximation 
needs no further comments, so we move on to the calculation of the total electromagnetic 
torque acting on the droplet within the PSA approximation. 



A. Total external electromagnetic torque 

We assume that a focused light beam passes through a liquid crystal droplet and that the 
average molecular director inside the droplet is oriented perpendicular to the beam axis. We 
choose a reference system in which the z axis coincides with the beam axis and the average 
molecular director lies in the xy plane. Note that, even if initially the average director of the 
droplet will not necessarily lie in the xy plane, the electromagnetic torque itself will force 
it there, in order to align the director to the optical electric field. So, at steady state, our 
assumption will be always verified. 



In order to calculate the total electromagnetic torque we must use either Eq. ( 14 ) or Eq. 



(15), with a = em and with expressions (jTj) and ^ of the torque density and stress tensor, 
respectively. The main difficulty is that the field to be used in the integrals is the total one, 
including both the external input field and the diffracted or scattered one. Neglecting the 
latter will give a vanishing result. So we need to calculate the propagation of light in the 
birefringent droplet. The first approximation introduced here for this calculation consists 
of simply replacing the droplet with a uniform planar slab of nematic liquid crystal having 
the same molecular director as the average one in the droplet. Therefore, all diffraction 
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effects are neglected and the only optical effects left to be considered are the changes of 
polarization due to birefringence (and dichroism) and eventually the attenuation due to 
absorption. Moreover, any distortion of the director configuration induced by light itself 
is assumed to be negligible, due to elastic interactions. As a second approximation, we 
treat the input light as a monochromatic plane wave propagating along the slab normal 
z. These two approximations combined are here named "planar symmetry approximation" 
(PSA). We stress that, despite its common usage [22l |23l [2l], PSA is very rough for the 
typical experimental geometry of strongly focused light beams and rather small droplets. 
Therefore, we can only anticipate semi-quantitative accuracy of its predictions. For example, 
in a strongly focused beam a large fraction of optical energy is actually associated with waves 
propagating obliquely, at a large angle with respect to 2;, which will see a much reduced 
birefringence with respect to the PSA plane wave. At any rate, all the model inaccuracies 
associated with the PSA approximation will not be very different for pure and dye-doped 
droplets. 

The slab thickness is taken equal to the droplet diameter d = 2R. The liquid crystal 
birefringence is denoted as An = ng — Uo, where no = Re (v^) and Ue = Re {\^e± + Ea) are 
the ordinary and extraordinary refractive indices, respectively. The absorption coefficient is 
denoted as ao (we neglect the dichroism for simplicity). The input light beam properties 
are the total power Pq, angular frequency 00 = 2ttc/X, vacuum wavenumber k = 27c/X, and a 
polarization assumed to be elliptical with its major axis parallel to the x axis and a degree 
of ellipticity fixed by the reduced Stokes parameter S3 or equivalently the ellipsometry angle 
X (in a complex representation of the input plane wave, their definition is S3 = sin(2x) = 
2lm{E^E;)/{\E^\^ + \Ey\^)). 

The calculation of the output wave fields emerging at the end of the slab is lengthy 



but straightforward, so we skip it. Inserting the input and output fields in Eq. (15) and 
integrating (the integration surface dV will be given by the two planes delimiting the slab 
and corresponding to the input and output fields; moreover, it is necessary to start the 
calculation with a finite wave and then take the plane- wave limit only after having performed 
a first integration by parts [25]), we obtain the following final expression of the external 
electromagnetic torque [22] : 

^em _ ^ r ^ h _ cos(A0)" 

00 ^ 
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1 - s?^ e~^"o^ sin(A0) sin , (18) 



where 6 is the angle between the director n and the x axis within the xy plane and A0 
2kRAn is the total birefringence retardation phase. 



It is interesting to note that the two main terms appearing in Eq. (18) tend to induce 

conflicting dynamics. The first, maximized for a circularly polarized input light (53 = ±1) 

and independent of the director orientation, tends to induce a constant rotation around the 

z axis, in a direction fixed by the sign of S3. The second term, instead, maximized for a 

linearly polarized input light (33 = 0) and dependent on the director orientation, tends to 

align the average molecular director of the droplet either parallel or perpendicular to the 

X axis, i.e., the major axis of the input polarization ellipse, depending on the birefringence 

retardation A0. For small values of S3 the latter term dominates and there is always an 

equilibrium angle 6 at which M|™ vanishes. If, instead, the polarization ellipticity |s3| is 

larger than a certain threshold s^t = sin(2xf) such that the former term dominates for any 

value of 6, the torque M|™ cannot vanish and the droplet must keep rotating (although not 

uniformly, unless S3 = ±1). The threshold ellipticity is defined by the following equation: 

S3t . ,^ . e-^°»^sin(A0) 

tan(2x0 = o^„ff„„ (19) 



I- s% 1 - e-2"ofl cos(A0) ■ 



It will be useful to consider also the mathematical limit of Eq. (18) for ao 00. This 
corresponds to the case in which the torque contribution of the light emerging from the 
output plane of the slab vanishes completely, as it occurs for very large absorption. However, 



we note that the same mathematical result is also obtained by taking the average of Eq. (18) 
over a wide range of birefringence retardations A(/), so that oscillating terms are canceled 
out. Such an average may occur as a result of two factors neglected in our PSA model: 
(i) oblique propagation of strongly focused light in the droplet (in our opinion this is the 
strongest effect), leading to reduced An and hence A^; (ii) propagation of light off droplet 
center, leading to an optical path length that is shorter than 2R. Whatever the actual cause, 
in this limit the electromagnetic torque reduces to the simple expression 

= Ml^ = — , (20) 

CO 

corresponding to the total flux of "spin" angular momentum associated with the input light 
only. 

Let us now turn to the droplet dynamics. 
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B. Droplet dynamics in the rigid body approximation (RBA) 



The RBA approximation can be justified by the fact that the typical viscosity (Leshe's) 
coefficients of the nematic liquid crystals (a typical value is 71 ~ 100 cP) are much larger 
than the water viscosity (r/ ^ 1 cP at room temperature). So any internal shear or relative 
rotation of the director with respect to the droplet fluid will be much slower than the overall 
droplet rotation with respect to the surrounding water. Moreover, it is possible that rigid- 
body behavior of the droplet (in particular in the steady-state dynamical regimes) is further 
enforced by the elastic interactions in combination with anchoring conditions (this second 
effect is especially plausible in the case of imperfect sphericity of the droplets). 

Within the RBA, the fluid velocity in the droplet is given by 

v(t) = ri(t)xr, (21) 

where Q{t) is the droplet angular velocity. Moreover, the molecular director is taken to 
rotate, everywhere in the droplet, at the same angular velocity as the fluid, i.e., it satisfies 
the following equation 

^ = f2(t)xn(t). (22) 

Since the orientation of the average director of the droplet in the xy plane is given by the 
angle 6 introduced in the previous section, one also has fi^ = dO/dt. 

As with all rigid bodies, all one needs in order to determine the droplet rotational dy- 



namics is Eq. (12) for the angular momentum rate of change, as applied to the entire droplet 



volume V = V^. By introducing Eqs. (21) and (22) in Eq. (13), the total angular momentum 
of the droplet can be rewritten as 

L = /O, (23) 

where / is the total moment of inertia, given by / = /y^ p(x^ + ip)dV = 2mi?^/5, and we 
have neglected /„. 

The total external torques M" acting on the droplet can be computed more conveniently 



using their surface integral expression, Eq. ( 15 ). As already discussed in the previous Section, 



the surface integral can be actually evaluated on a surface that lies just outside the liquid 
crystal droplet boundary, i.e., within water, thereby making the calculation much simpler. 
Let us now consider each of the five contributions a = hp, el, vis, em, ph. 
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First, owing to the spherical shape of the droplet (within the SDA), the pressure torque 
M^P will vanish identically, as it can be readily verified by a direct calculation. Since in 
water there are no elastic stresses, the elastic torque M*^' will also vanish identically. 

The viscous term does not vanish and it can be easily evaluated by solving the Navier- 
Stokes equations in water with assigned velocity on the droplet boundary as given by Eq. 



(21) in the laminar fiow limit (and neglecting the effect of the electromagnetic stresses in 
water). The result of such calculation is the well known Stokes formula for the rotational 
viscous torque acting on a rotating sphere in a viscous fluid: 

M™ = -6r]Vdn, (24) 

where rj is the water viscosity coefficient. 

The electromagnetic torque M''™ does not vanish. Within the PSA approximation dis- 



cussed above, its z component will be given by Eq. (18), while its x and y components will 
vanish. 

Finally, as already discussed in the previous Section, the photoinduced torque Mp*^ should 
also vanish based on our theory, as in water there should be no photoinduced stresses. 
However, as discussed above, we cannot exclude that a higher-order theory might predict 
a nonvanishing M^^ associated to interfacial effects (a specific possibility for such an effect 
would be, for example, a photoinduced discontinuity of the fiow velocity at the boundary 
between liquid crystal and water). Therefore, we should consider this possibility in our 
analysis. 



Neglecting inertial terms, equation ( 12 ) in the RBA model is then reduced to the following 
torque balance: 

M™ + M^"' + = 0, (25) 



where the first two torques are given by Eq. (24) and (18), while the expression of Mp*^ is 



unknown. Equation (25) is actually a first-order differential equation in the rotational angle 

e{t). 

Let us assume initially that Mp*^ = 0. As discussed in the previous section, the steady- 



state solution of Eq. (25) depends on the value of the polarization ellipticity S3 with respect 



to the threshold value s^t given in Eq. (19). For [ssl < s^t the solution is static, i.e 



6'(t) = Oo is constant, while for |s3| > sst the solution is dynamical and corresponds to a 



14 



generally nonuniform rotation of the droplet around the z axis. In the circular polarization 
limit the rotation becomes uniform. 

By a simple integration, it is possible to determine the overall rotation frequency / of 
the droplet, which takes the following expression: 

,2 



/ = /o Re I (l - e-2"°^cos 

-(l-s^)e-4"o«sin2^^]i|^ (26) 

where 

/o = Po/(167r2a;r7i?3) (27) 



(note that Eq. (26) includes also the stationary solutions / = 0, for |s3| < s^t)- The highest 
frequency is obviously reached for S3 = ±1, i.e., for circular polarization of the input light. 
Note also that if we take the ao ^ cxd limit (which may be actually due to all the factors 



discussed above and neglected in PSA), we obtain simply / = /o instead of Eq. (26). 

Let us now consider what should happen instead for Mp'^ 7^ 0. As we said, we do 
not know the actual expression of a nonvanishing Mp'^, as this should result from some 
unknown higher-order term in the constitutive equations. However, since the photoinduced 
torque density r^*^ acting on the molecular director is proportional to the electromagnetic 
one T*^™, it is reasonable to expect that also this photoinduced torque Mp'^ is proportional 
to M*^™. The ratio of the photoinduced to electromagnetic torque density is C/^a, a number 
which is of the order of several hundreds. The corresponding ratio of total external torques 
is therefore limited by this value, although it could be smaller. 

In the case of circularly polarized input light, the effect of the photoinduced external 
torque would be that of inducing a dye-enhanced droplet rotation, as revealed by a higher 
rotation frequency achieved for the same input light power, or a smaller light power needed 
to obtain the same rotation frequency when compared with the undoped case or to what 
happens when light falling outside the dye absorption band is used. In the experimental 
section we will specifically search for such effects. 



C. Uniform director approximation (UDA) 



According to the RBA model presented in the previous subsection, the rotation speed of 
a droplet is independent of all photoinduced effects, unless higher-order interfacial effects 
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should be found to be significant. This conclusion relies strongly on the assumption that 
the droplet rotates effectively like a rigid body. In this subsection we analyze theoretically a 
situation in which the director is not constrained to rotate together with the fluid flow. Such 
a situation may occur with appropriate boundary conditions for the director orientation and 
fluid motion at the droplet interface. We show that in this case the velocity field and the 
director rotation are influenced essentially by photoinduced effects. 

In order to keep the model manageable analytically, we assume here that the spatial 
distribution of the director in the droplet is always approximately uniform (this assumption 
was unnecessary in the RBA model). Moreover, we will still consider the fluid motion to 
coincide with that of a rigid body. 



The dynamics of v and n flelds will therefore still be taken to be given by Eqs. (21 ) and 



(22), but the two fi's entering these equations will be different, in general. Let us then label 
fly and fin the angular velocities of the fluid and director rotations, respectively. 

We can flnd two dynamical equations for f2„ and f2„ starting from Eqs. ([T]) and (|2]), 
respectively, and following a procedure similar to that used when finding Eq. ( [l2| . First, we 
integrate both sides of Eq. ([T| over the whole droplet volume Vd- Second, we multiply both 
sides of Eq. ^ vectorially by r and integrate over the whole volume Vd- Owing to Eq. (16), 
we can then write the two resulting equations in the following form: 

inVd^ = 

dt 

= M*°* - T'°\ (28) 



where the total external torque M*°* = J2a is still defined by Eq. (15) and it is therefore 
identical to that used in the RBA model, and we have also introduced the total internal 
torque T*°* = J2a exchanged between director and velocity degrees of freedom. More 
precisely, for each kind of interaction a = hp, el, vis, em, ph, the torque T" is simply defined 
as the volume integral of the torque density r" over the whole droplet. 



Equations (28) highlight the coupling between the droplet fiuid rotation and the director 
dynamics inside the droplet. One recovers the RBA model when the internal torque T'°* is 
very rigid and acts as a constraint that locks the director to the fiuid. 

Since ~ with excellent accuracy, from the first of Eqs. (28) one finds T*°* ^ 0. 



Then the second of Eqs. (28) becomes identical to that of the RBA. This result, however. 



should not be taken as saying that within the UDA the droplet always behaves exactly as 
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in the RBA, as the external torques M", and in particular the electromagnetic one, can be 
affected by the director orientation in the droplet, which will not be the same as in RBA. 
In particular, as we will see, the photoinduced effects will be present even in the first-order 
theory, in which they only appear via the internal torque T^^. 

Let us now calculate the internal torques T" using the explicit constitutive dependence of 
the corresponding director torque densities r". The hydrostatic pressure contribution van- 
ishes identically. For a perfectly uniform director the elastic torque density t'^^ and therefore 
the total torque T*'^ also vanish identically. However it is not obvious that one can truly 
neglect the elastic torque density, as it is just this torque that keeps the director approxi- 
mately uniform. If all other torque densities are always uniform then the director will remain 
uniform "spontaneously" , and the elastic torque density will indeed vanish. Whenever the 
other torque densities are nonuniform (in our case, for example, the electromagnetic torque 
density is not very uniform in large droplets), their effect will give rise to a nonvanishing 
elastic torque density which balances them, in order to constrain the director to remain 
approximately uniform. At any rate, we assume here that the total elastic torque, obtained 
by integrating this elastic torque density over the whole droplet, remains always negligible. 
Further work will be needed to assess the validity of this assumption. 

Let us now consider the other interactions. The viscous torque under our assumptions is 
easily calculated and is given by 

T^" = -7iVd(fl„-ri.), (29) 



where 71 is the orientational viscosity coefficient [Ti 

As for the electromagnetic term, under our hypothesis of uniform director (UDA) and 
spherical droplet (SDA) of constant density, it can be shown that the following identity holds 
with high accuracy: 

The proof of this identity is reported in the Appendix. Thereby, the calculation of T°™ = 



M^™ can be based on the PSA approximation, and its explicit expression is given by Eq. ( 18 ). 

Finally, the photoinduced term (in first-order theory, the only one we consider in this sec- 
tion) is simply related to the electromagnetic one when the latter is caused by an optical wave 
having a wavelength within the dye absorption band. In this case, since r^^^ = (C/ea)T'^™, 
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the photoinduced internal torque will be exactly given by 

= ^M*^'". (31) 

However, it must be kept in mind that this strict relationship between the photoinduced and 
the electromagnetic torques is only valid when the eletromagnetic field is that of an optical 
wave absorbed by the dye. It is instead possible to separately tune the electromagnetic 
and photoinduced torques by adding a second wave whose wavelength is out of the dye 
absorption band, or alternatively by adding static electromagnetic fields. In these cases, 



Eq. (31) will only apply to the contribution of the wave having a frequency within the dye 



absorption band. 



By using these results, Eqs. (28) can be rewritten in the following dynamical equations, 
which define our UDA model of the droplet dynamics: 

InVj^ = -7iV,(f2„ - ri,) + M^- + ^ 0, 

= -6Vrfr/n, + 7iV,(f2„ - n,) - (32) 

where in the second equation we have cancelled the two electromagnetic terms M''™ and 
— T'^™ as they are exactly opposite to each other. 



Equations (32) show explicitly the internal nature of the photoinduced effects (according 
to first-order theory), with the "action" and "reaction" torques appearing respectively in the 
former and latter equation (or vice versa). In contrast, the electromagnetic torque is external 
(the "reaction" term acts on the electromagnetic field) and is applied on the director only, 
in the first place. It is only via the viscous internal coupling between the director and the 
fiuid motion (and possibly also the elastic one, which we neglected) that the electromagnetic 
torque finally drives also the droplet fiuid rotation. 



Let us now study the dynamics predicted by Eqs. (32) in a couple of interesting examples. 

Let us assume first that a single circularly polarized laser beam travelling along the z 
axis drives the droplet rotation. The director will spontaneously orient in the xy plane. 
Then, the electromagnetic torque becomes independent of the director orientation within 



the plane, so that its rotation may become uniform at steady state. In this case, Eqs. (32) 
are readily solved and give 

" V. + 71 J ' 
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^v. = ^: (33) 

where M^™ can be approximately calculated using Eq. (18). These results show that the 
droplet fluid and the director do not rotate at the same rate. For C = (no photoinduced 
effect), the difference in angular velocity is very small as the viscous coefficient 71 is about 
two orders of magnitude larger than rj. In contrast, in the presence of photoinduced effects, 
the ratio C/^a can be larger than a few hundreds [31 [26], so that a significant difference in 
the two angular velocities should become possible. 

As a second example, let us assume that there are two linearly polarized optical waves 
travelling along z and driving the droplet dynamics: one is assumed to be polarized along 
X and to have a wavelength falling outside the dye absorption band; the second is instead 
assumed to be polarized at an angle \E' with the x axis and to have a wavelength that is 
within the dye absorption band. Having two waves with only one being absorbed by the 
dye, the electromagnetic and photoinduced torques can be adjusted independently to each 
other. 

Being linearly polarized, the two waves will induce torques that tend to align the droplet 
director parallel or perpendicular to the polarization direction of the respective wave. So 
we may assume that at steady state the director will acquire a fixed orientation along some 
intermediate direction 6 with < 6' < \E'. Therefore, at steady state one may assume 
fin = and the electromagnetic and photoinduced torques to be constant in time. By 



solving Eqs. (32) with these assumptions, one finds that the equilibrium orientation of the 



director is actually fixed by the following balance: 

+ — = 0, (34) 
71 + 6// ^ ^ ^ 

while the droplet fluid rotates constantly with the following angular velocity: 

= ^ = — ' , . (35) 

This result is rather counterintuitive: a continuous droplet rotation is induced by two linearly 
polarized optical waves having different planes of polarization. This could not occur without 
the presence of photoinduced effects, i.e., based only on the electromagnetic torque (it must 
not be forgotten that the two waves have different wavelengths, so their superposition is 
not coherent). We see, then, that the photoinduced effects associated with the dye can give 
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rise to highly nontrivial effects in a hquid crystal droplet, as long as the director can be 
decoupled from the fluid motion. 

IV. EXPERIMENT 

We prepared emulsions of pure and dye-doped liquid crystal (LC) in bidistilled water. 
By properly choosing the relative quantities of liquid crystal and water, we could obtain 
relatively stable emulsions containing many liquid crystal droplets having micrometric size, 
most of them in the 1-20 fim diameter range. Most experiments were performed using the 
commercial liquid crystal mixture E63 (see Ref. [27j for its composition), provided by Merck, 
Darmstadt, Germany. This material is convenient for the wide temperature range of its ne- 
matic phase (from -30°C to 82°C). A few experiments were performed using the liquid crystal 
4-cyano-4'-pentyl-biphenyl (5CB). The dye used in doped materials was the 1,8-dihydroxy 
4,5-diamino 2,7-diisopentyl anthraquinone (HK271, provided by Nematel, Mainz, Germany), 
known as one of the most effective dyes in the photoinduced effects. We prepared dye-liquid 
crystal solutions at a concentration of 2% in weight, leading to torque enhancement ratios 
(/sa of several hundreds [3l [26]. 

For performing the optical rotation experiments, a small volume of emulsion was placed 
on a microscope glass slide and covered with a thin (thickness of about 0.18 mm) glass 
coverslip, thus forming a thin cell (open on the sides). A 100 x oil-immersion microscope 
objective (Carl Zeiss, NA=1.25) was used both to image the droplets on a CCD camera 
(using lamp illumination from below the cell, in combination with other optics) and to focus 
the input laser beams that were used to trap the droplets and induce their rotation, as 
discussed in the following (see Fig. [T]). The objective was attached, via an oil thin film, to 
the glass coverslip. When needed, the microscopic imaging could be made in the crossed 
or nearly-crossed polarizers geometry, in order to visualize the droplet birefringence. CCD 
images and movies (acquisition rate of 50 frame/s) could also be recorded on a PC, for 
subsequent frame-by-frame analysis. 

The setup for the optical trapping and manipulation is a dual-wavelength optical-tweezers 
apparatus (see Fig. [T]). Our setup differs from usual optical tweezers arrangement in that 
two laser beams, respectively generated by a diode-pumped solid state laser (A =785 nm, 
subsequently called "IR" beam) and a He-Ne laser (A=633 nm), are brought to a common 
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focus at the specimen plane of the microscope objective. Both beams were able to achieve 
trapping of droplets for diameters in the 2-12 /im range, as proved by a sudden stop (or 
confinement) of the droplet Brownian motion after trapping. The beam-waist (Airy disk) 
radius at focus is estimated to be about 0.4 fim, in the adopted objective-overfilling configu- 
ration. Actually, due to the thickness of the glass coverslip, our optical trap center (roughly 
corresponding to the laser beam waist) could not be located right in the middle of the cell 
but was close to the covershp. This fact might have led to some contacting of the droplets 
with the glass, particularly in the case of the largest droplets. In such cases the droplets 
may have experienced a somewhat stronger friction, and therefore our results on the rota- 
tional speed could be biased. However, any systematic effect will be identical for pure and 
dye-doped droplets and obviously independent of the light wavelength used for inducing the 
rotation. So, our comparison will remain valid. Moreover, the agreement between theory 
and experiments will show that these possible systematic effects are certainly small, at least 
for not too large droplets. 

The common focus of the two lasers allowed us to switch the controlling beam from one 
wavelength to the other (e.g., to switch the He-Ne beam on and the IR off, or vice versa) 
without changing the trapped liquid crystal droplet. Since the He-Ne wavelength is close to 
the maximum of the dye absorption band, while the IR wavelength falls completely outside 
the absorption band, this corresponds to adding or removing the photoinduced effects, so 
as to allow for a direct comparison of the behavior obtained in the two cases with the same 
droplet. 

Besides for trapping, the two beams were also used to induce the electromagnetic and pho- 
toinduced torques acting on the droplet and driving its reorientation or continuous rotation. 
By inserting along the beam (just before the microscope objective) a suitable birefringent 
wave plate, we could control the polarization of each beam (separately). In particular, we 
have used linear polarizations with an adjustable polarization plane and elliptical polariza- 
tions with an adjustable degree of ellipticity S3 = sin(2x). We could also obtain certain 
specific polarization combinations of both beams together, such as one linear and the other 
circular or elliptical. By using a suitably dispersive optically-active plate (home-built us- 
ing a chiral-doped randomly oriented liquid crystal cell) we could also obtain two linearly 
polarized beams with two different polarization planes, forming an angle of ^ ~ 40°. 

Micron-sized LC droplets are known to show either one of two possible director spatial 
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distribution [201 Ell 121]: axial (or bipolar), already discussed in the previous section, and 
radial, which has full spherical symmetry and a hedgehog defect at the droplet center. The 
images of axial and radial droplets under a polarizing microscope can be almost identical. 
However, we could identify axial droplets by checking that their image changed if the mi- 
croscope polarizers were rotated, or by looking directly at their dynamical behavior under 
the laser beam, as only axial droplets could be readily set in rotation if illuminated with 
circular polarizations or reoriented using linear polarizations. Radial droplets could not be 
rotated at all (because all the optical torques vanish, due to radial symmetry, if absorption 
is neglected). About 80-90% of the droplets of our emulsions were found to be axial In 
the following we will only refer to them. 

Using linearly polarized light, we could easily fix the average director orientation of our 
trapped axial droplets. By looking at the microscopic image pattern under crossed polariz- 
ers [291 EQ], we verified that the director orientation was indeed parallel to the polarization 
plane. By slowly rotating the polarization by means of a half-wave plate, the droplet align- 
ment followed the polarization. These experiments could be done both with pure liquid 
crystal droplets and with dye-doped ones, the latter both with IR and He-Ne beams, with- 
out much difference. 

Using elliptically polarized beams and axial droplets, we could perform the analogue of 
the classic Beth's experiment [31], i.e., set the trapped droplets in continuous rotation by 
the transfer of angular momentum with light. Depending on the droplet size and on the 
laser beam power and ellipticity, the rotation could range from very slow (periods of sev- 
eral seconds) to very fast (down to the millisecond range). We used two different methods 
for measuring the spinning of droplets by analyzing the frames of a droplet rotation movie. 
First, using the crossed-polarizers geometry we could image the birefringence rotation, corre- 
sponding to the director motion, as a periodic pattern modulation in time. For the smallest 
droplets, we could observe a significant modulation of light intensity during the rotation 
even without the analyzer, presumably due to the anisotropic scattering cross-section of the 
droplets. Second, we could also measure the droplet fluid rotation (i.e., regardless of director 
rotation) when a smaller satellite droplet (or another small object) happened to get trapped 
close to the rotating droplet and was thereby dragged around, as in the case shown in Fig. |2| 
The rotation speed measured in these two ways (when the satellite was sufficiently small 
and close to the rotating droplet) were always found to be the same. 
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In order to check if the photoinduced effects give rise to a rotation speed enhancement, on 
each droplet we measured the rotation frequency induced by the He-Ne and the IR beams 
as a function of input beam polarization ellipticity. We repeated this for many different 
droplets having a range of diameters. We also measured the rotation frequencies induced 
in pure (undoped) liquid crystal droplets, although of course in this case the comparison 
could not be done with exactly the same droplet sizes. In all these measurements, the two 
laser beams were adjusted for having a roughly equal angular momentum flux, as given by 



Eq. (20). In particular, the light power measured after the microscope objective was about 
4.1 mW for the He-Ne beam and 2.8 mW for the IR in most data shown. These values 
correspond, for a given polarization ellipticity S3, to the same angular momentum flux to 
within 15% (M|™ = (1.3±0.1) x 10^^^ Nm for the case of circular polarization). However we 
also investigated the power dependence of the rotation frequency in some droplets. 

An example of the measured rotational frequencies of a given droplet versus light po- 
larization ellipticity (given by the angle x) is shown in Fig. |3} We actually measured this 
dependence for many other droplets of different sizes, made of both pure and dye-doped 



LC. Each of these measurements was then fitted by means of Eq. (26). In these fits the 
radius R was fixed at the value determined by analyzing the droplet microscopic picture. 
The absorbance at 633 nm was instead measured separately on a bulk sample (we obtained 
Oq = (1.0 ± 0.2) X 10^ cm~^) and then kept fixed to this value in all fits. The constants 
/o and A0 were adjusted for best fit. /o was found to be roughly consistent (to within a 



factor two in most cases) with its predicted value, as given by Eq. (27) (using the known or 
measured values of the laser power and frequency, water viscosity and droplet radius) . From 
the best-fit value of A0, we could also estimate the birefringence An, which was almost 
always found to fall in the range 0.11-0.13 (average 0.12), with a few droplets giving a value 
of 0.10 and 0.14. These values are inconsistent with the known refractive index difference of 
the bulk material {An = 0.2273 at 589 nm, Merck data sheet). We ascribe this discrepancy 
to the strong approximations associated with the PSA model (in particular to neglecting the 
contribution of obliquely propagating waves in the focused beam), as a strong depression of 
the optical anisotropy due to confinement effects is not plausible for micrometric droplets. 
It is apparent from Fig. |3] that the rotational behavior of a dye-doped LC droplet under 
IR or He-Ne laser beams is not identical. However, the difference is well explained by the 



different wavelength and absorbance in Eq. (26), while no particular rotation enhancement 
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is seen in the He-Ne case, where photoinduced dye effects should tal<;e place. 

In Fig. |4] we show the rotation frequency observed for several pure LC (open symbols) 
and dye-doped LC (filled symbols) droplets of different sizes using circularly polarized light. 
Panels (a) and (b) refer to rotations induced by the IR and He-Ne beam, respectively. 
The insets show the corresponding (linear) dependence on laser power for a fixed droplet 
(incidentally, this linear behavior supports our assumption that the laser light induces no 
significant distortion of the director configuration in the droplet, in the power range used in 
our experiments). 

Again, a first conclusion that can be immediately drawn from this figure is that no 
significant rotation speed enhancement takes place in the dye-doped LC case with the He- 
Ne beam as compared both to the IR beam case and to the pure LC case. Should the 
photoinduced torque be external, one would have expected a rotation speed enhancement 
by a factor of the order of C/^a, i-e., of several hundreds. 



In Fig. |4j the dashed lines are the predictions of Eq. (26) after adjusting the laser power 
P for best-fit to data. The birefringence was kept fixed to the average value An = 0.12 
obtained from the measurements discussed above (and confirmed by the threshold ellipticity 
data, as discussed below), but increasing its value led to worse fits (in particular. An = 0.23 
leads to very bad fits). From Fig. |4| it is seen that the agreement between data and theory 
is reasonable, although there is a statistically significant discrepancy. In particular, the data 



do not show at all the oscillations predicted by Eq. (26) (of course, small residual oscillations 
might be hidden in the noise). Moreover, for both lasers the best-fit values of the power 
were found to be about a factor two smaller than the actual measured values (assuming 
a water viscosity rj = 1 cP, corresponding to the room temperature of 20° C). The solid 



lines in Fig. |4] correspond instead to the simpler theory f{R) = /o(-R), as given by Eq. (27), 
obtained in the limit ao oo. In this case, no adjustable parameter was used, i.e., the values 
of the laser power and of the water viscosity are both fixed to the known values [rj = 1 cP). 
Nevertheless, it is seen that the agreement is even better and that the same theory (taking 
into account the difference in light power and frequency) explains all the data, i.e., both the 
dye-doped case with He-Ne beam, when there is significant absorption, and the IR beam 
or pure LC cases, when there is no absorption (incidentally, this agreement shows that any 
systematic effect due to the droplet closeness to the glass wall is essentially negligible, except 
perhaps for the largest droplets). This better agreement of the simplified theory obtained 
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for tto — ^ oo clearly cannot be truly ascribed to absorption (negligible in the IR and pure 
droplets case). It is instead the likely result of light diffraction, oblique propagation and 
other effects neglected in the simple PSA calculation, which may all contribute to averaging 
out the oscillations due to the outgoing light, as discussed previously. This view is also 
confirmed by the fact that previous works on transparent liquid crystal droplets reported 
similar observations [231 El] • 

Let us now consider the behavior of the threshold ellipticity Xt for droplet rotation versus 
the droplet radius R. For each droplet and wavelength, Xt was obtained from the best fits 
performed on the measured rotation frequency versus light ellipticity (such as that shown in 



Fig. In particular, xt is entirely determined by the best-fit value of A0, via Eq. (19). The 



resulting data are shown in Fig. |5| together with the predictions obtained using Eq. (19) 
combined with A0 = 2kRAn and adjusting the birefringence An for best fit. It is seen 
that the theory agrees reasonably well with the experiment. The best fit is for An ^ 0.12 
for both laser wavelengths, confirming the average value obtained before. Moreover, for the 
He-Ne case the effect of light absorption is clearly seen in the data. 



It must be noted that Eq. (19) is based on the PSA model including the effect of the 
light emerging from the droplet. Indeed, taking the ao ^ oo limit leads to = for all 
droplet sizes, which is inconsistent with our data. We do not know exactly why the PSA- 
neglected effects do not affect the threshold ellipticity as much as the rotation frequency 
(apart from the depression of An), but this is what we actually find and what was also 
reported in previous works [231 [21J. At any rate, we stress that all inaccuracies of PSA 
theory do not affect our comparison between pure and dye-doped droplets or between IR 
and He-Ne beams. 

Finally, we tried to find evidences of a decoupling between the droplet fluid and the 
average director degrees of freedom. In particular, we searched for the effects predicted by 
the UDA model to occur when the He-Ne and IR beams are simultaneously impinging on the 
droplet, both linearly polarized with the two polarization planes forming an angle of about 
40°. In order to visualize the droplet fluid rotation independently of the director, we imaged 
the dynamics of droplets having small satellites (such as smaller droplets or other particles 
in the suspension), so that the droplet surface flow would be highlighted by the revolution 
of these dragged objects. However, despite many efforts, we could not find any sign of a 
steady-state decoupling between the director and the fluid motion in the droplets (except, 
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perhaps, for some complex transient effects) and, tlierefore, of a corresponding continuous 
fluid rotation of the droplet. This probably indicates that elastic effects associated with the 
nonuniform director distribution and a not perfectly spherical shape of the droplet (effects 
neglected in the UDA model), concur to keeping the orientational and flow degrees of freedom 
locked to each other. 

V. CONCLUSIONS 

In summary, we have studied for the first time the dynamical behavior of droplets of dye- 
doped nematic liquid crystal trapped in water by a strongly focussed laser beam and set in 
rotation by the optical torques generated by the same beam. We searched for evidences of a 
dye-induced enhancement of the droplet rotation speed associated with photoinduced effects, 
but no enhancement was found. This null result is in accordance with the leading first-order 
theory of these photoinduced effects [H E] and directly proves that the dye-enhanced optical 
torque is not associated with an exchange of angular momentum with light or other external 
degrees of freedom, but that it must instead be associated with a fully internal exchange 
of angular momentum between the molecular orientation and fluid flow degrees of freedom. 
The latter exchange is correctly described only by assuming the existence of a dye-induced 
stress tensor acting on the fluid flow in the illuminated region, as first proposed in Ref. 
[5]. This photoinduced stress tensor embodies the internal "reaction" to the photoinduced 
torque acting on the molecular director. Our null result is therefore also a direct proof of 
the actual existence of this stress tensor. 

The interplay between fluid flow and director dynamics within the droplet is also pre- 
dicted by our models to give rise to very peculiar effects, such as a continuous droplet fluid 
rotation induced by two linearly polarized optical waves having different polarization planes. 
However, the occurrence of these phenomena requires to unlock the constraint existing be- 
tween the droplet fluid and the molecular director. This constraint is partly due to viscous 
forces, but it may also be due to elastic interactions combined with anchoring forces. Ex- 
perimentally, we could not reach a situation in which this constraint was effectively broken, 
so we could not demonstrate the predicted effect. 
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APPENDIX: PROOF OF THE ELECTROMAGNETIC TORQUE IDENTITY 



To prove the validity of identity (30), we start from Eq. (14) as applied to the case 



a = em, which can be rewritten as follows: 

M'r-Tr = h= 1 6.,.xAC- (A.i) 



V 



We must therefore prove that the integral li vanishes identically within the SDA and UDA 
approximations. 

Let us consider the first term in Eq. ([8| of the electromagnetic stress tensor. This is an 



isotropic pressure-like term, i.e., of the form p6hk- For such a term, the integral (A.I) can 
be recast into a surface integral as shown in the following: 



/ eijhXjdkpShkdV = / eijhXjdhpdV 

Jv J dV 



IV JdV 

Within the SDA, the surface element normal dAh and the position vector Xj will be parallel 
to each other, and therefore their cross product (as expressed by their product times the 
antisymmetric Levi-Civita tensor) vanishes identically. 

Let us now consider the second term in Eq. ([s]) of the stress tensor. Omitting the time- 
average for brevity, we must consider the following integral 



r(2) 



X e.,,xMDuE,)dV = I e.,,x,DAE,dV 
= eijhXjDkdhEkdV - eijhXjDkekhi^dV 
= / ^ijhXjSkiEidhEkdV — I ^ijhX jDkekhi-^dV, 

where we have used the two Maxwell equations dkDk = and dkEh — dhEk = —tkhidBi/dt 
and introduced the dielectric tensor to express D in terms of E. The dielectric tensor is 
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always symmetric for a permutation of its indices. Moreover, within UDA and constant 
density approximations, the dielectric tensor is also uniform within the droplet. Exploiting 
these two properties, we obtain 

dB, 



r(2) 



1 /" f uB 

,dh{eMEiEk)dV- jDkekhi-r^dV 



where the first integral vanishes identically (within SDA) as it has taken a pressure-like form. 

The third stress-tensor term, in the magnetic field, can be treated analogously to the 
second, obtaining the following expression: 

dD, 



j(3) 



Summing up the three contributions, the whole integral = + jf^^ + if^ is reduced 
to the following: 

d 



eijhXj-^{€khiDiBk)dV 

dV. (A. 2) 



rx|(DxB) 

This expression (which is a part of the electromagnetic torque due the so-called Abraham 
force [13], which is completed if one considers also the effect of the neglected electromagnetic 
force mentioned in [16]) corresponds to a true residual electromagnetic torque which does 
not vanish in general. However, after optical-cycle averaging, its order of magnitude can be 
estimated as follows: 

hc^\B\\E\vJ—)c^Tr(—), (A.3) 



r / V r / 

where r is a characteristic time of the droplet dynamics and At ^ R/c is the droplet crossing 
time of light. In our case At/r ^ 10"^^, making this contribution totally negligible with 



respect to Tf™ and therefore the identity (30) valid to high accuracy. 
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FIG. 1: Experimental setup. Legend: P polarizer, WP half- wave plate, DM dichroic mirror, L lens, 
CCD - ccd camera, BE beam expander for objective back apertme overfilling, X either retardation 
wave-plate (for the appropriate laser wavelength) or optically active plate (see text) for rotating 
the linear polarizations of the two wavelengths by a different amount, depending on the specific 
experiment. 
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FIG. 2: Sequence of photograms showing the hght-induced rotation of a LC droplet in time, as 
highhghted by the revolution of a dragged small object. Scale-bar: l^um. 
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FIG. 3: Droplet rotation frequency / versus input light ellipticity angle X) for a dye-doped LC 
droplet having a radius of 2.4 ^um, using the IR (a) or He-Ne (b) laser beams, with a power of 
2.8 mW and 4.1 mW, respectively. The dots are the measured values and the solid line is the 



theoretical fit based on Eq. ( 26 1 . The threshold ellipticity xt is obtained from the fit from the two 
symmetrical points at which the solid line crosses zero. 
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droplet radius R (|iin) 

FIG. 4: Frequency / of the droplet rotation induced by circularly polarized light. The main panels 
show the dependence on droplet radius for a fixed laser power (power on the sample: 4.1 mW for 
the He-Ne and 2.8 mW for the IR beams), while the insets show the dependence on laser power for 
a fixed droplet radius of 1.8 /um. Panel (a) refers to rotations induced by the IR laser, panel (b) 
to the He-Ne case. Filled (open) circles refer to droplets made of dye-doped (pure) liquid crystal. 
Solid and dashed lines are the theoretical predictions obtained as explained in the text. 
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FIG. 5: Threshold ellipticity for droplet rotation xt as a function of droplet radius. Data points 
are actually obtained from the fits described in the text and in the caption of Fig. |3j The solid line 



is from Eq. (19). Panel (a) refers to the IR case, panel (b) to the He-Ne. All data are for droplet 
of dye-doped liquid crystal. 
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